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00 ' An exactly soluble model describing a spin 5=1 antiferromagnetic chain doped with mobile 

S — 1/2 carriers is constructed. In its continuum limit the undoped state is described by three 
ON 1 gapless Majorana fermions composing the SU(2) triplet. Doping adds to this a scalar charge field 

and a singlet Majorana fermion with different velocity. We argue that this mode survives when the 
Haldane gap is added. 

< 

It is well known that mobile holes introduced into antiferromagnetic Mott insulators by means of doping cause 
frustration which seriously affects the magnetic properties of the system. In this paper we study a doped spin-1 chain 
which in the presence of mobile holes interpolates between S = 1 and S = 1/2 states. According to Refs. jj]U this 
situation is realized in the carrier doped Haldane system Y2- a; Ca a ;BaNi05 For x = the interactions between 
the spin-1 Ni 2+ ions are well described by a Heisenberg model, upon doping mixing of the 5=1/2 holes on the 
Oxygen sites leads to a low energy doublet state in an effective one-band Hamiltonian which can move in the S = 1 
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background. 



The model we study is an integrable model whose Hamiltonian closely resembles the Hamiltonian introduced in [Q : 

L 



ri=l 



H%™ = -(l-5 SilSj )V iJ (S t .S j ) . 

> 

ly-j • Here Sf = Si(Si + 1) with Si = 1 or 1/2 and Vij permutes the spins on sites i and j. Comparing (jl|) to the effective 
■^j- ' one-band Hamiltonian for the doped Ni oxide given in g], 
t-H ! 

q ■ = S Si s 3 ,iJSi ■ Sj+i - Vij (Si ■ Sj + - j (2) 

we find several differences. First, the spin exchange between 5=1 sites contains biquadratic terms giving the spin-1 
Takhtajan-Babujian chain for hole concentration x = 0. Therefore in the undoped limit the spectrum is gapless and 
one may think that the most spectacular feature of the doped S = 1 chain - filling the gap with holes, is lost in our 
model. We shall see later however, that it is possible to reintroduce a gap in the continuous limit where we have 
"-^J ■ a field theoretical description of the model. Second, due to a hidden supersymmetry of the model (see below) the 
ratio between the exchange integral and the hopping is fixed in (^). We do not expect this to be a serious restriction 
on physical properties of the model. Third, hopping processes such as |lj ij) — > | [i lj) are allowed. Finally the 
integrable model contains an antiferromagnetic exchange interaction between nearest neighbour hole states which is 
responsible for the absence of a ferromagnetic phase of (||) for sufficiently large hole concentration. Such an interaction 
has, however, been considered in Ref. || to improve agreement with the experiments at x = 1. 

In the following we will use the exact solution of ([!]) to deduce an effective field theory for the low-energy sector of 
the system which will enable us to study possible relevant perturbations. The low-energy limit of the undoped system 
(x — 0) is well known to be a realization of an SU(2) level k = 2 Wess-Zumino-Novikov-Witten (WZNW) model 
with central charge c = 3/2. This model is equivalent to a model of three massless Majorana fermions composing an 
SU(2) triplet (qj. A completely filled (x = 1) band corresponds to the S = 1/2 chain which in the low-energy limit 
is equivalent to the SU(2)i WZNW model. For finite doping we find, as may be expected, one free bosonic mode in 
the charge sector. The spin sector however, turns out to be rather unusual containing a direct sum of c = 3/2 and 
c = 1/2 models with different velocities. Thus doping generates the fourth Majorana fermion - a feature observed in 
various models related to two-channel Kondo physics |(||7j]. 

The model (Q) is constructed from solutions to the Yang Baxter (YBE) equation invariant under the action of the 
graded Lie algebra gl(2\l) in the 'atypical' representation [S} + M. These representations contain two multiplets of 



1 




spin S and (S — 1/2) with respect to the 5'?7(2)-subalgebra of gl(2\l). A well known example for a lattice model 
obtained from the [l/2] + -representation in this approach is the supersymmetric t-J model (see e.g. ||). Denoting 
the generators of the S'C/(2)-subalgebra by S a , the [/(l)-charge operator by B and the remaining fermionic generators 
of gl(2\l) by V^, (which create and annihilate the holes, see ||) the local ^-operator of the Quantum Inverse 
Scattering Method reads 

p + 2iB 

C(p) oc | i^2V+ p + i(B + S z ) -iS+ \ . (3) 
-iy/2V~ 

The spectrum of the transfer matrix for a vertex model with L of these weights per row is obtained by means of 
the algebraic Bethe Ansatz px|| . Starting from the fully polarized state of spin-S multiplets on each site we consider 
states with Nh holes (generating sites with spin S — 1/2) and magnetization M z = LS — \Nh — This leads to an 
auxiliary eigenvalue problem for an inhomogenous graded six- vertex model on a lattice of n — Nh + iVj_ sites which is 
solved by means of a second Bethe Ansatz through addition of states with fermionic grading (holes) to its eigenstate 
with n — iVj_ (for details see e.g. Q). As a result the spectrum of the supersymmetric vertex model is parametrized 
by the roots of the following set of Bethe Ansatz equations (BAE) 

Aj + ID \ _ T-r Aj — Afc + l -pr Aj — v a — 5 

j^l,...,N h + N l (4) 
1= x 7 ' a = l,...,N h . 

For S — 1/2, where the model becomes the supersymmetric t-J model, and Eqs. (Q) are Sutherland's form of the BAE 
p| . The grading of the underlying algebra allows to construct two more sets of (equivalent) BAE [^2|,^). Eqs. ([|) are 
most convenient for our analysis, however. 

To construct a local Hamiltonian for S — 1 we project the graded tensor product £ a (fi+ ^)®£b(/x— |) in two copies 
V a , Vb of the matrix space of (^) onto the [1]+ representation in V a ® Vb- The result £[ 1 ] + (/x) of this fusion [^3 14 



defines a new vertex model with eigenstates again parametrized by the roots of (|4|). The logarithmic derivative of its 
transfer matrix at \i — is the Hamiltonian (Q) with eigenvalues 

E{{\ 3 },{v a })-HM Z -fiN h 

Nh+Nl f 2 \ Nh ( 1 \ 

= H ( H -xrn)-^{» + \ E )- LE < 5 » 

fe=l V fc 1 / Q = l v / 

(we have added an external magnetic H field and a (hole) chemical potential jl = fj, + 1 to the Hamiltonian) . 

In the thermodynamic limit the solutions of the BAE (^) are real hole rapidities v, whose density we denote as p{y) 
and complex A-strings with densities cr n (A) (n = 1, 2, . . .). The equations for the densities are 

5 n ,2s(a;) = c7„(ie) + C nm * a m (x) - 5 n ,is * p{x) 

(6) 

a 2 * s(x) = p(x) + [1 + a 2 ] * p(x) + s * a\(x) 

where all densities marked by tilde correspond to holes, f * g{x) denotes a convolution, s(x) = 1/(2 cosh ttx), a n (x) = 
{2n/n)/(Ax 2 + n 2 ) and 

C nm (x) = 5 nm 5(x) - (8 n +\ tm + S n -i,m) s(x) 
Writing the energy (||) in terms of densities of A-holes and p we obtain 

E/L = -2 - J dx [27r(a 2 * s) + p] p{x) 

+ / dx2TTs(x)a2(x) — lim Hn / dxa n (x) . (7) 

J rwoo J 
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By minimization the free energy we obtain the thermodynamic Bethe ansatz (TBA) equations for the energies 
e n = Tln((7 n /(T„) of A-strings and k = Tln(/5/p) for the hole rapidities 

e n {x) = Ts * ln[l + e'-'W^p + e ^+^)/ T ] 

-27rS ni2 s(x) - 5 nA Ts * ln[l + e~ K ^' T ] , (8) 

subject to the condition lim, i ^ oc (e n /ri) = H and 

-[27ra 2 * s(x) + fi] - Ts * ln[l + e ei(a;)/T ] 

= k(x) + TR * ln[l + c - k{x)/t ] (9) 

where R = a 2 * (1 + a%) ■ In terms of these functions the free energy reads 

—T { dx{a 2 * s)(x) ln[l + e - Ki - x ^ T ] . (10) 



For temperatures T <C H the energies e„> 2 can be eliminated from (||), giving the following equation for e 2 

e 2 (x)-Ti?*ln[l + c e2(2:)/T ] 

= -2?rs(a;) + H/2 + Ts * ln[l + e ei{x),T ] (11) 

Together with the n = 1 equation from (^) and (^) this equation determines the low temperature phase diagram of 
the system (Fig. ^). For sufficiently large negative /i Eq. (^) implies n > corresponding to vanishing hole density. 
The remaining TBA coincide with those for the integrable S = 1 magnet fll5]]l6|l . Choosing fi > H/2 we find n < 
and ei < which implies a hole density x = 1. The resulting TBA for e„> 2 are those of the 5 = 1/2 Heisenberg 
chain. Similar considerations for intermediate /i yield the complete phase diagram and lead us to identify k and e± 
as the energy of the charge and spin modes associated with the mobile carriers while e 2 is the magnetic mode of the 
background spins. 

We now concentrate on the case of finite doping in a vanishing magnetic field which corresponds to chemical 
potentials n — 4 < fi < 0. In this regime we have k < for < Q. For temperatures smaller than the Fermi 
energy of holes T <C —fi one can replace k in Eq. (|^) by its zero temperature value ko(x) and the free terms by their 
asymptotics. As a result we get 

e n ( x ) = Ts * ln[l + e^-W 7 ) [i + e ^+i(x)/T^ 

-^Tr^e-^l - 2nA6 nA e- 7r ^ (12) 

where A = — (2-k)^ 1 J®q dye 7ry . To study the specific thermodynamic properties of the present model we have to 
separate the contributions to the free energy stemming from the charge-sector from those due to the e n . Considering 
low temperatures again the leading contributions to k come from the vicinity of the Fermi wave vectors ±Q. In this 
region one can safely neglect contributions to Eq. (^]) from e\ and rewrite it as 

-[27ra 2 * s(x) +n]-TR* ln[l + e ~ |K(a:)l/T ] 

rQ 

= k(x) - / dyR(x - y)n(y) (13) 



where Q is determined by the condition k(±Q) = 0. Using the procedure introduced by Takahashi jT7[ , we can rewrite 
the free energy ( |To| ) as F/L = Eq /L + fi + / 2 where E is the ground state energy and 

fi = -T J dx P n{x) ln[l + e -\ K ^ x )\/ T ] « - w T 2 /6v , (14) 
-T [ dx(s * p )(x) ln[l + e ei ^/ T ] . (15) 
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Here po(x) is the density of hole rapidities at T = 



p Q (x)- dyR(x - y)po(y) = a 2 * s(x) 

J-Q 

and v is the velocity of the charge mode k 

' ■ (16) 



2np G (Q) dx 



In the low-T limit T <C —fi where this form is valid one may replace s{x) and s * po{x) in ( p"5| ) by their asymptotics. 
Now we have TBA equations in such form that the spin sector is manifestly decoupled from the charge one. The 
free energy of the latter is given by Eq. (|l4|) representing a scalar bosonic mode, while Eqs. ( |l2|) and (|15| ) describe 
the thermodynamics of the spin sector. Such TBA equations (i.e. with two driving terms) arise also in other systems 
combining different spins fl8f| and in the two-channel Kondo model with a channel anisotropy [ fl9|| . As shown in 
Ref. JtJ such two-channel Kondo model can be written in terms of four Majorana fermions. These considerations lead 
us to hypothesize that the effective low-energy theory for the Hamiltonian ([!]) at small doping is: 

H cB = J dx [H + H] (17) 

. 3 . 

"H = - 2 X! V 2Xad x Xa ~ 2 VlX ° dxX0 + -9X0X1X2X3 
a=l 

and a similar expression Ti({x\) for the left movers. The fields \0i Xa are Majorana fermions and the marginal 
coupling g between two sectors in (^) is the only one possible by symmetry (we have omitted marginally irrelevant 
interactions between sectors with different chirality). A magnetic field couples to the term i J (flxoXa + ^abcXbXc) 
quadratic in the fermions with some constant (3. 

The parameters of the effective theory have to be extracted from the TBA at small doping. This corresponds to 
A -C 1 and implies that the two terms in ( |l5| ) are dominated by contributions from the region 7rx ~ ln(27r/T) (where 
| £2! ~ T) and by ttx ~ \n(2irA/T), respectively. This separation of scales allows to obtain the small A behaviour of 
( |l5| ) at low temperatures 

h = -*Jl(l-ZA lnA )^Jl(z + zA lnA 



For A — this is the free energy of a single Majorana fermion (xo i n (P-7[)) with velocity v\ and that of the SU{2)2 
WZNW model with velocity v 2 which can be expressed in terms of massless triplet of Majorana fermions ||. The 
ratio of the velocities is given by 

1 I%^yn {y) 

V1/V2 = -— — . (18) 

The coupling constant g needs to be chosen such that it produces the ^4-dependence of j 2 - Similarly, the parameter 
(3 can be determined from the T = magnetic susceptibility. At small A -C 1 we get 

X=-[ —\nA+. 

IT \V2 Vl Z7T 

giving (3 2 = A/2tt. In Fig. ^ some of these quantities are given as a function of the hole concentration x. 

The simplest S'f7(2)-invariant perturbation of Jl7j ) is a mass term i J dx (m Y]^— 1 XaXa + MxoXo^j f° r the particles 
in the spin sector (recall that the model of three massive Majorana fermions is a good continuous approximation for 
the S = 1 Heisenberg antiferromagnet [^0)). Choosing M = we get a situation with a Haldane gap in the triplet 
sector and gapless singlet excitations. This may be a good description for the magnetic states inside the Haldane gap 
observed in the doped Ni oxide M . 

To summarize, we have introduced an integrable model describing a magnetic system which interpolates between 
the integrable S = 1 Takhtajan-Babujian and the S — 1/2 Heisenberg chain. For finite doping we find up to three 
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massless modes determining the low temperature thermodynamics of the system. In the SU (2) invariant case (H = 0) 
these are a bosonic charge mode and a direct sum of a SU (2)2 WZNW model and one Majorana fermion with different 
velocities in the spin sector. Further analysis of the TBA is necessary to gain more insight in the nature of the phase 
diagram of this and the related S > 1 systems. 

This work has been supported by the Deutsche Forschungsgemeinschaft under Grant No. Fr 737/2-2. 
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FIG. 1. Low temperature phase diagram of the doped chain in a magnetic field H 3> T: the saturated ferromagnetic phase 
is labelled FM. At intermediate fields the low energy properties of the system are determined by gapless bosonic modes with 
energies k (A, Bi), ei (A, B2, FM) and £2 (A, Bi,2, C). The H = phase is discussed in the main text. 
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